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Abstract 

The inhomogeneous distribution of matter in the non-linear regime of 
galaxies, clusters of galaxies and voids is described by an exact, spheri- 
cally symmetric inhomogeneous solution of Einstein's gravitational field 
equations, corresponding to an under-dense void. The solution becomes 
the homogeneous and isotropic Einstein-de Sitter solution for a red shift 
z > 10 — 20, which describes the matter dominated CMB data with small 
inhomogeneities. A spatial volume averaging of physical quantities is in- 
troduced and the averaged time evolution expansion parameter 9 can give 
rise in the late-time universe to a volume averaged deceleration param- 
eter (q) that is negative for a positive matter density. This allows for a 
region of accelerated expansion which does not require a cosmological con- 
stant. A negative deceleration parameter can be derived by this volume 
averaging procedure from the Lemaitre-Tolman-Bondi open void solution, 
which describes the late-time non-linear regime associated with galaxies 
and under-dense voids and solves the "coincidence" problem. 



email: john.moffat@perimeterinstitute.ca 

1 Introduction 



In recent articles, we investigated a cosmology in which a spherically sym- 
metric inhomogeneous enhancement is embedded in an asymptotic FLRW uni- 
verse 1 The inhomogeneous enhancement is described by an exact inhomo- 
geneous solution of Einstein's field equations. The inhomogeneous distribution 
of matter at late times can possibly allow for a negative deceleration parameter 
q, depending on the nature of the exact inhomogeneous solution of Einstein's 
field equations. This will allow for the possibility of explaining the acceleration 
of the universe, without a cosmological constant. The model also leads to an 

1 For a more complete list of references, see refs. HI 121 
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axis pointing towards the center of a spherically symmetric inhomogeneous en- 
hancement with dipole, quadrupole and octopole moments aligned with the axis. 
The distribution of CMB temperature fluctuations can be unevenly distributed 
in the northern and southern hemispheres f . 

We shall describe the inhomogeneous late-time universe by an exact spher- 
ically symmetric solution of Einstein's field equations. The solution becomes 
the homogeneous and isotropic Einstein-de Sitter matter dominated solution 
for z > 20, so that it can describe the WMAP data near the surface of last scat- 
tering with small temperature perturbations 5T/T ~ 3 . However, for 
z < 20 the solution is inhomogeneous and for luminosity distances correspond- 
ing toz~ 0.2 — 1.4, includes the Type SNIa supernova measurements We 
investigate the expansion geometry of the spherical distribution of inhomoge- 
neous matter and under-dense voids to see whether it allows for the possibility 
of a negative q and an accelerating expansion of the universe, without negative 
pressure dark energy or a cosmological constant. 

A generic relativistic solution of a three-dimensional inhomogeneous universe 
is presently unknown. The highly symmetric Lemaitre-Tolman-Bondi (LTB) 7 
8 9 solution that we shall consider in the following is inhomogeneous in only 
one spatial dimension. Even though the assumption of spherical symmetry is 
unrealistic, the model is useful in illustrating important local features of physical 
quantities. 

For our inhomogeneous late-time universe model, we are required to carry 
out a volume averaging of physical scalar quantities, such as the expansion 
parameter 8. We find for such an averaging process that for a zero cosmo- 
logical constant and for irrotational matter the average deceleration parameter 
(q) can be negative and describe an accelerating local region of the universe, 
corresponding to our observed Hubble radius, without invoking a cosmological 
constant. 

2 Inhomogeneous Friedmann Equations 

For the sake of notational clarity, we write the FLRW line element 



where k = +1,0,-1 for a closed, flat and open universe, respectively, and 
dil 2 = d0 2 + s'm8 2 d<fi 2 . The general, spherically symmetric inhomogeneous line 
element is given by [3 El El HDl IH1 C21 C31 IU : 




(1) 



ds 2 = dt 2 - X 2 {r,t)dr 2 - R 2 (r,t)dfl 2 . 



(2) 



The energy-momentum tensor T^ v takes the barytropic form 



T» u = {p + P )u»u, 



(3) 
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where w' 1 = dx^/ds and, in general, the density p — p(r,t) and the pressure 
P = p( r i t) depend on both r and t. We have for comoving coordinates u° — 
1, u % = 0,(i= 1, 2, 3) and g^ v u^u v = 1. 
The Einstein gravitational equations are 

G„ v + k 9ilv = SwGT^, (4) 

where G M „ = R^ — ^g^TZ, 1Z = g^R^u and A is the cosmological constant. 
Solving the Goi = equation for the metric we find that 

, \ R'(r,t) 

x{r ' t)= fW' (5) 

where R' — dR/dr and f(r) is an arbitrary function of r. 
We obtain the two generalized Friedmann equations pQ: 

i? li? 2 11., ... 1-R'i? 1//' 4ttG, „. 1. 

i^+S^ + S^^- 7 ) "3^i? + 3i^? = - — (P + 3P)+ 3 A ' (?) 
where R = dR/dt. 



3 Late-Time Matter Dominated Universe 

The late-time matter dominated universe will be pictured as a large-scale inho- 
mogeneous enhancement that is described by an exact inhomogeneous spheri- 
cally symmetric solution of Einstein's held equations. The inhomogeneous en- 
hancement is embedded in a matter dominated universe that approaches asymp- 
totically an Einstein-de Sitter universe as t — > oo. An observer will be off-center 
from the origin of coordinates of the spherical inhomogeneous enhancement. 
The local inhomogeneous solution corresponds to an open, hyperbolic solution 
and describes an under-dense void. Surveys such as the 2-degree Field Galaxy 
Redshift Survey and the Sloan Digital Sky Survey show large volumes of rela- 
tively empty space, or voids, in the distribution of galaxies [H]. 

For the matter dominated Lemaitre-Tolman-Bondi (LTB) [71 |S1 E] model 
with zero pressure p = and zero cosmological constant A = 0, the Einstein 
held equations demand that R(r, t) satisfies 

2R(r, t)R 2 (r, t) + 2R(r, i)(l - / 2 (r)) = F(r), (8) 

with F being an arbitrary function of r of class C 2 . 
The proper density of matter can be expressed as 

F' 

9 = WirGR'R 2 ■ ^ 
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We can solve (0 to obtain 
where 

H^ = H 2 ± + 2H ± H r , (11) 
is an effective Hubble parameter and 

* = f, *x-f. (12) 
We have for the critical density 

8ttG Pc = H 2 S . (13) 

There are three possibilities for the curvature of spacetime: 1) f 2 > 1 open 
(fi - 1 < 0), 2) / 2 = 1 flat (fi - 1 = 0), / 2 < 1 closed (fi - 1 > 0). 

We describe the inhomogeneous density regime by a hyperbolic, open solu- 
tion f 2 (r) > 1 with p = A = 0, corresponding to an under-dense void. We 
choose 

f(r) = \/l + r 2 (14) 

and the metric reduces to 

ds 2 = dt 2 _ R*ir,t) dT 2 _ R 2^ ^2 (15) 

A parametric solution is given by 

i?(r,i) = lF(r)(f 2 (r) - l)- 1 (cosh u(r,t) - 1), (16) 

t + [3(r) = ±(f 2 (r)-l)- 3 / 2 (smhu(r,t)~u(r,t)) 7 (17) 



**<Wr,t) = ( ^ ) ( ) ( J - (18) 



-2^(r)\/ / 2 (r)-l ' 
F 2 (r)^ i?'(r,i) ; Vsinh 4 i U (r,t), 

Here, /3(r) is an arbitrary function of r of class C 2 JU| . 

The function /3(r) can be specified in terms of a density on some spacelike 
hypersurface t = to- The metric and density are singular for 

u = or R'(r,t)=0, (19) 

where u — refers to the hypersurface t + (3(r) ~ and the singular surface 
for R'(r,t) = is more complicated. Our pressureless model requires that the 
singular surface t(r) = £(r) describes the surface on which the universe becomes 
matter dominated (in the FLRW model this occurs at z ~ 10 4 ). The hyperbolic 
FLRW model is obtained by choosing the conditions 

/ 2 (r) = l + r 2 , F(r) = 4sr 3 , (3{r) = 0, (20) 
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where s is a positive constant. This leads to the metric solution 

R = sr(coshu- 1), (21) 
t = s(smhu-u), (22) 

-(— 

4s 2 Uinh 6 



iirGp = ( ~T~TE~r~ ) ' (23) 



We shall change the radial coordinate to 

r* =r- r SLSl (24) 

where r = tsls is the location of the surface of last scattering corresponding to 
r* = at a red shift z ~ 1100. We postulate that (3{r*) approaches a maximum 
constant value at the approximate end of galaxy and cluster formation, Zf ~ 
1 — 5, where z = Zf denotes the final value of the red shift. We assume that 
P(r*) — > and f(r*) — > 1 as r* — > 0. A model for j3 is of the form: 

ft r * ) = a (bT?Y' (25) 

where a, b and a are positive constants, (3(0) = and (3(r*) — > a for r* ^> &. 
The surface r* = for the metric is a singular surface. 
The metric for f(r*) — ► 1 becomes 

ds 2 = di 2 - (t + /3(r*)) 4 / 3 (F 2 (r*, i)dr* 2 + r* W), (26) 

where 

ir c-"- 1 + i^y < 27 > 

and 

1 

6nG(t + fi(7*)) 2 Y(r*,t)' 

We obtain for j3(r* ) — > the homogeneous and isotropic solution corresponding 
to a spatially flat Einstein-de Sitter universe with the metric 

ds 2 = dt 2 - a 2 (t){dr* 2 + r* 2 dn 2 ), (29) 

and a(t) oc t 2 / 3 . This solution is compatible as a background spacetime with 
the approximately homogeneous and isotropic WMAP data for 20 < z < 10 3 
and Sp/p - 10" 5 0. 

The Einstein-de Sitter solution must be supplemented by non-zero pressure 
contributions for z ~ 10 3 in the linear regime at the surface of last scattering in 
order to take into account the radiation density pr. Moreover, a more general 
solution is required that has a non-vanishing cosmological constant A, so that 
we have Q = Qm + fijj + £l\ = 1 in agreement with the spatially flat result 
£1 = 1.02±0.02, obtained from the determination of the first acoustic peak in the 
WMAP power spectrum 3 . However, in the following section, we shall replace 
f^A by a contribution corresponding to an acceleration parameter q obtained 
from our inhomogeneous void solution with A = 0. 
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4 The Inhomogeneous Late-Time Density and 
Deceleration Parameters 



We obtain from ijfjj for f(r) = y/l + r 2 and p — (we replace for convenience 
r* by r): 

t2 , r; ' _ , \ (30) 



//: +2H r H x -^-^ = 8nGp + A. 



Dividing this equation by H\ gives 



n 



inGp A 



2r 



(31) 



Let us expand R{r, t) in a Taylor series 



R(r,t) = R[r,t -(t -t)} = R(r,t ) 



l-(t -t) 



f?(r,t ) , 1 



iJ(r,*o) 2 



+ ^(to-*)^ (Mo) 



R(r,t ) 



R(r,t ) 



1 - (to - t)Bb± - 2 (*o ~ *)«(»•, *o)Hoj 



(32) 



where t denotes the present epoch and H j_ = R(r,t )/R(r,t ). Moreover, we 
have 

R(r,t )R(r,t ) 



(33) 



R 2 {r, t Q ) 

For the matter dominated era we set p — in Eqs.® and Q and substitute 
for R from J7J to obtain 

1 4irp Q {r,t ) 



l(r, to) = r + 00-2 / 

3 3H$ ± (r,t ) 



A 



lF 0r (r,t ) 1 r 2 



3RR'H$ ± 3H$ x (r,t ) 3H ±(r,t ) 3R 2 H 2 A 



(34) 



where H 0r (r,t ) = R'(r,t )/R'{r,t ). 

If we set H ±(r,to) = Ho r (r,to) = If (to) where If (to) = a(to)/a(to) and 
i?(r, to) = a(to), then we obtain the global FLRW expression for the deceleration 
parameter for a spatially flat universe with f 2 (r) = 1: 



9o — 2^ 0M — ^ 0A ' 



(35) 



where from (|31|l f2oM + ^oa — 1 with CIqm — 8irpoM /3H 2 (to) and flo\ = 
A :!//-.:/.;,). 

The variance of q is given by the exact non-perturbative expression: 



var( 9 ) = (q 2 - (q) 2 } 



2\i/2 _ 



q, 



(36) 
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where 

— — f d 3 x./j(...) , 

(...) = J rj 1 ^' , 37 

denotes the ensemble average and 7 denotes the determinant of the 3-dimensional 
metric gij (i = 1, 2, 3). 

We see from (1341) that different observers located in different causally dis- 
connected parts of the sky will observe different values for the deceleration 
parameter q, depending upon their location and distance from the center of 
the spherically symmetric inhomogeneous enhancement. This can lead to one 
form of cosmic variance, because the spatial average of all the observed values 
of local physical quantities, including the deceleration parameter, q, will have 
an intrinsic uncertainty. 

Let us set the cosmological constant to zero, A = 2 . From (|31 \ and (|34|) . 
we obtain for t = tg: 



n = + f + 2r 2Hor = 1, (38) 

H QX R 2 Hq ± R'RHq ± H ± 



and 



1 4nGf> (r,t ) 1 r I H Qr (r,t ) 1 r 2 

q[r ' t0) 3 + 3flg ± (r,t ) 3RR'H* ± 3H 0± (r,t ) 3i? 2 ff 2 / ^ 



'0_L 

If in l|39|) we have 

2 

H 0r H 01 _ + T — > AirGpa + H 2 ± + ^ , (40) 

then the deceleration parameter q can be negative and cause the universe to ac- 
celerate without a cosmological constant or dark energy. If we solve for Hq t / Hoi_ 
from Ea. l|3^|) and substitute into i|39[l. we obtain 

^=5-K^r)' (41) 

which gives a negative q for 

r 2 /R 2 H^ > 3/5, (42) 

and we can satisfy the condition |J5SJ. We see that for r = corresponding to 
f 2 (r) = 1 and a spatially flat universe q(r,t ) = 1/2. For p = A = 0, Eqs.|0 
and (|40l) lead to R > as the condition for an accelerating expansion of the 
universe. 

The conditions l|38|) and H4U|) require that f 2 (r) > 1 for the inhomogeneous 
enhancement, corresponding to an under-dense void. If we choose instead to 
describe the inhomogeneity by the spatially flat solution f 2 (r) — 1, then we 
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We do not provide any solution to the cosmological constant problem, namely, why A = 0. 
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cannot simultaneously satisfy SI = 1 and q < 0. However, we will investigate 
in the next Section, whether the condition l|40|) can indeed be satisfied for the 
particular LTB model that we have adopted. We shall find that a more gen- 
eral inhomogeneous solution is needed to establish whether such a solution can 
explain the accelerating expansion of the universe without a negative pressure 
fluid or a cosmological constant. 



5 The Raychaudhuri Equation 

The Raychaudhuri equation is [T4*| 



dO 1 
— = -R uv V> i V v + 2oj 2 - 2a 2 - -6 2 , (43) 
as 3 



where we have adopted geodesic world lines, is a time-like vector and 



UJ 



2 - u^Up, > 0, o 2 — o^g^ > 0. (44) 



Here, o ^ v and u)^ denote the shear tensor and the vorticity tensor, respectively. 
From Einstein's field equations, we have for p = A = 0: 

^ = 2cj 2 - 2a 2 - -9 2 - iirGp. (45) 
as 3 

The deceleration parameter q is defined by 

(3dd/ds + 9 2 ) 6(cr 2 -uj 2 ) +4irGp 



e 2 e 2 



(46) 



If the vorticity uj — 0, then the deceleration parameter q has to be positive for 
positive p. 

For our synchronous comoving coordinates we have 



dO _ 
ds 

We have 

= (30 + e 2 ) _ 

q ~ e 2 

We may now conclude that if 



V^^O = V^d^O = 0. (47) 



3 |1 + 6 | + 6 2 - ZH 2 - m\ 



(48) 



3H 2 + 6H 2 ± > 3^-+6^- + 9 2 , (49) 

then q could be negative leading to a local acceleration of the universe. However, 
for our synchronous comoving frame, we have 

uj^ = V [V V M] = d v V^ - dfjYu = 0, (50) 
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so that the vorticity for this chosen gauge is zero. Therefore, we must deduce 
that the condition (|49|l cannot be satisfied for the positive energy condition 
p > and p > 0. Moreover, the two conditions (|40|l and (|42|l cannot be satisfied 
for the synchronous and comoving metric However, in the following section, 
we will show that we must carry out a spatial volume averaging of the expansion 
parameter 9 and its time evolution to arrive at a physically viable description 
of the local inhomogeneous late-time universe and an accelerating expansion of 
the universe. 



6 Spatially Averaged Cosmological Domains 

For our inhomogeneous model, it is necessary to perform a spatial average l|37[) 
of physical quantities, due to their observer, location dependence. Let us define 
a more specific domain averaged expression for a scalar quantity |15l 1161 117| : 

(*(f,t))u - rj- / d 3 x^(x,t), (5i) 

Vv Jd 

where 

Vd = [ d 3 x^f (52) 

is the volume of the simply-connected domain, D, in a hypersurface. We can 
define an effective scale-factor for our spatially averaged spherically symmetric 
inhomogeneous enhancement (void): 

(R(r,t)) D =(V$°) 1/3 , (53) 



where ViD is the initial volume. 

The volume averaging of the scalar \1/ does not commute with its time evo- 
lution dsi mini: 

{Mr,t))D-dt{*{r,t)) D = (9(r,t)) D {9(r,t)) D - (*{r,t)6{r,t)) D . (54) 
We can derive for the averaged (9) d the equation 

(9> D = 3j^=3(ffx)D. (55) 
{H)d 

Setting ^ = 9 and substituting for the Raychaudhuri equation for an irrotational 
matter dominated late-time model with lj^ = 0, d6/ds = d9/dt and A = 0, we 
obtain the expression 

3 77^ = -^G(p)d + (P)n, (56) 

where (P)d is a function of (9)d, the average Hubble expansion parameters 
(H±) d and (H r ) d and the average shear (a) jj. 
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For inhomogeneous cosmology, the smoothing due to averaging of the Ein- 
stein field equations does not commute with the time evolution of the non-linear 
field equations. This leads to extra contributions in the effective, averaged Ein- 
stein field equations, which do not satisfy the usual energy conditions even 
though they are satisfied by the original energy-momentum tensor. It is the 
lack of commutativity of the time evolution of the expansion of the universe 
in a local patch inside our Hubble horizon that circumvents the no-go theorem 
based on the local Raychaudhuri equation 

If we can satisfy the condition 

(P) D > 4nG(p) D , (57) 

then the averaged deceleration parameter (q)D can be negative corresponding 
to the acceleration of a local patch of the universe. If we substitute the value 
of (R)d from the volume averaged equation for A = 0, then the exact 
LTB solution for an irrotational inhomogeneous, under-dense void solution can 
describe the local acceleration of the late-time universe and fulfill the constraint 
equation 1)57(1. The spatial domain average of q in Ea. (|41|l gives 

This yields the averaged condition for an accelerating universe: 

> 3/5. (59) 



This spatial domain averaged condition for our LTB void solution can be satis- 
fied in our synchronous comoving gauge with zero vorticity for (p)d > 0. 



7 Conclusions 

We have modelled the late-time inhomogeneous and non- linear regime for z < 10 
by an exact spherically symmetric, inhomogeneous solution of Einstein's field 
equations. The solution approaches an Einstein-de Sitter matter dominated 
solution for z > 10, which describes a homogeneous and isotropic FLRW back- 
ground spacetime with small inhomogeneities Sp/p ~ 10~ 5 near the surface of 
last scattering in agreement with the WMAP data A volume averaging 
of physical quantities is necessitated for our inhomogeneous universe, and due 
to the lack of commutativity of the time evolution of the average of the scalar 
expansion 9, we find that it is possible to have a negative averaged deceleration 
parameter (q) for a positive density of matter p and a zero vorticity ui = 0. Ob- 
servational bounds on the magnitude of the quantity (Pd) in the condition l|57|l 
must be obtained to establish whether a sufficiently negative q can be produced 
by the expanding inhomogeneous enhancement or void to allow an explanation 
of the accelerating universe without a negative pressure fluid or a cosmological 
constant. 
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The inhomogeneous solution is exact, so that we do not have to be concerned 
with the failure of perturbative backreaction calculations in the non-linear late- 
time universe regime. However, we have assumed a high degree of symmetry 
for our inhomogeneous solution, since it only allows for a one-dimensional in- 
homogeneity in the radial direction, although the angular azimuthal expansion 
described by H± does play a significant role in determining the cosmological 
solution. It would be interesting to investigate inhomogeneous cosmological 
solutions that possess less symmetry than the LTB models |13|. 
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